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ABSTRACT 


One topic of current interest in thermoacoustic research is an annular prime 
mover [Lin et al. J. Acoust. Soc. Am, 100, 2846 (1996)]. The starting point for 
this research is an investigation of a constricted annular resonator. A literature 
search of the field resulted in surprisingly few references. The results of analytic, 
numerical, and experimental investigations are presented. Introducing a constriction 
into an annular resonator splits each longitudinal duct mode into two modes, one of a 
higher frequency with a pressure antinode at the constriction and one at a lower 
frequency with a velocity antinode near the constriction. The lower mode is more 
sensitive to changes in the length and porosity of the constriction than the higher mode. 
Overall agreement between measured and predicted mode shapes and resonance 
frequencies is very good. It was found that it is necessary to include end 
corrections at the constriction to get accurate agreement between measured and 


predicted results. 
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L INTRODUCTION 


The purpose of this thesis is to investigate acoustic resonances in a constricted, 
annular resonator. The motivation for this project is based on a related problem in the 
area of thermoacoustis - the annular prime mover. The intent is that the measurements 
reported in this thesis will serve as preliminary validation of models being developed for 
the prime mover project. 

The azimuthal modes of a uniform cross section annulus are degenerate. The 
introduction of a constriction splits the degeneracy resulting in two modes having different 
resonance frequencies and distributions of acoustic pressure (mode shape). The goal of 
this thesis 1s to both measure and predict the frequency splitting and mode shapes. 

A theoretical discussion will be covered in the Chapter II. Chapter II is a 
description of the experimental apparatus including the annular acoustic resonator and 
peripheral equipment. In addition, a detailed discussion of the data collection procedure 
will be covered. Chapter IV is an analysis of the data obtained and comparison with 
predictions. The thesis is summarized and important conclusions are addressed in Chapter 
V. 





Il. THEORY 


A. INTRODUCTION 


The purpose of this chapter is to outline the theoretical analysis of a driven, 
constricted annular resonator. The discussion begins with the solution of the wave 
equation in and the effective radius of a uniform cross section annular resonator. This 
solution is developed in cylindrical coordinates. However, because the purpose of 
this thesis is to investigate frequencies below those required to excite cross modes, we 
adopt a simplified plane model for the majority of analysis. 

The next step is to introduce a constriction. It is shown that the constriction 
results in frequency splitting of what were degenerated, orthogonal modes in the uniform 
cross section resonator. Itis also shown that the degree to which a particular mode is 
excited depends upon the location of the driver relative to the constriction. Because of 
dissipation, the modes overlap in frequency space and, in general, two modes are 
simultaneously excited. When this situation occurs the resonance frequencies of the 
individual modes cannot be determined from peaks in the frequency response. 
However, pole-zeros analysis of the frequency response can yield accurate values for the 
resonance frequencies. The chapter concludes with a discussion of pole-zero analysis of 


two coupled modes. 


B. SOLUTION FOR A UNIFORM ANNULAR RESONATOR 


The analysis begins with the general solution of the acoustic wave equation 


in cylindrical coordinates 
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2n = 2.1 
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Assuming harmonic time dependence, the wave equation becomes the Helmholtz 


equation 
Oo 2 
V’P=-k’P where k* = (2) . (2.3) 


The general solution of Eq. (2.3) is assumed to be expressed as the product of three 


terms, each a function of only one variable. 


<{°)_ })x 


Figure 2-1. Geometry for the annular resonator. 


P(r,0,z) = R(r)O(@)Z(z) (2.4) 
Substituting Eq. (2.4) into (2.3) and dividing by ROZ , gives 




















ilegl. Oct ( ch) ,~t1 270. ,i28% = a OS) 

Rr ér\ @ér r? @ 80? Z éz? a 
Rearranging to isolate the z- dependent terms gives 

Li, 2k 1 1 6°0 ; le ae 

SS SS rs 2.6 

Rr or\’ dr)" +? © 80° Z és 
Following the standard separation of variables solution technique , 

ae 
let —_— = —-k?. oe 
° Z woz : ey 
The general solution to this equation is 
ZL=A, cos(k,z) +B, sin(k,z) (2.8) 


2 


Equation (2.6) now becomes after multiplying through by rr’ and 


rearranging 








Oo OR 7 
i2/, )+(e?-e2)?--242 : . (2.9) 





R @r Or ; C6? 
1 6’ 
Let ———— = -m’. . 
S ® oe? m (2.10) 
The general solution to this equation is conveniently expressed in the form 
© (8) = cos(m6+ ¢) (2.11) 
Equation (2.9) is now reduced to 
aon OR ) oo ee 
R Or rap + (k -r° —-m“)=0, (2712) 
where Kea ko Kk (2.13) 
The general solution to this equation is 
R(r)=A,J,(k,r)+B_Y (kr), (2.14) 


where J, (k,r) and Y. (k,r) are the cylindrical Bessel and Neumann functions. The 
Neumann function must be retained because the omgin is excluded in our annular 
resonator. 

Ail boundaries are assumed to be rigid, so the normal component of acoustic 
velocity must be zero at the boundaries. Applying rigid boundary conditions at z= 0 
and z= L, implies that B, =0 and 


A: | Sows, 20. (2.15) 
The function ©(@) must be periodical, in other words, ©(@) = ©(6 + 2z). This 
condition is satisfied if m is an integer, m=0,1,2 ----. 
Application of rigid boundary condition at r=aand r=b implies that 
A J m(k,a) +B _Y n(k,a) = 0 
A J m(k,b) +BY m(k,b) = 0. (2.16) 


Expressed in matrix form these equations become 


J(k,a) Y.({k,a)\(A,) (0 
He) eese)@ em 
The values of k,a that satisfy this equation are those for which the determinant of the 
first matrix is zero , 1.e. 
Joh @ Wire ioe b) Jka Dy ealheane=aOe (2.18) 
Let the values of k,a that satisfy Eq. (2.18) be called @,. With this definition, the 
first of Eqs. (2.16) becomes 


ces ak — 
and Eq. (2.14) becomes 
J. 
R(r) =A, [J,(4,") - eee (2.20) 


For the purposes of this thesis, we are interested in frequencies low compared 
to those required to excite cross modes. To eliminate pressure vanations in the z 
direction, we require that /=0 in Eq. (2.15). In this case Z(z)=1, and k? = k?’. 
We are also interested in modes for which the circumference of the resonator is 
approximately one wavelength long, i.e. m= 1 in Eq. (2.11). Under these conditions, 


the solution to the Helmholtz equation is 


P(7,O0)=A [J,(4,r) - Fea Y (ik, )]008(6 +o). (2.21) 


The dimensions for our particular resonator are a= 10.00 cm and b = 15.30 cm. The 
lowest value of a,, is then a, =0.796. Recall that for /=0, 


20 


he oe 
| e 





a a 
)*. Also, k, = a . Combining these expressions for k, 


implies that 


276 ea 


A a 


Let us define the effective circumference of the resonator L,, = mA, or for m= 1 


ée 





(2.22) 


L., = A . Combining this definition with the Eq. (2.22) and rearranging yields 


27a 
L eg = = (2.23) 
10 





™ . Using Eq. (2.23) , we 





if 
Alternatively, we can define an effective radius R,, = 5 





have RK, = . For @) = 0.796 , R,, = 12.56 cm. 
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Cc. PLANE WAVE THEORY 
The parameter R,, is defined in such a way that 


= —— = =—— = — 2.24 
Rey 27s ~ 27 oa) 
The angular coordinate 0 can be related to the arc length x measured at R,, through 


x 
G= KX, (2725) 
R egy 





With this definition , the @ - dependence of the acoustic pressure becomes 
cos(kx + ¢). (2.26) 
This function can be expressed in terms of complex exponentials as 
Gostkx +d) = 2fe 7) 4 e Tiet)) | (2.27) 
Substituting this result into Eq. (2.21) gives 
RitNOm b—mi@aniie ome" e fee ss’ 4| (2.28) 
The interpretation of Eq. (2.28) is that the pressure is represented by two counter 
propagating waves that behave essentially like plane waves. This analogy will be 


used to analyze the nature of the sound field in the annular resonator. 


The plane-wave-like solution will be valid for frequencies below those required 
to excite the first cross mode. The frequencies corresponding to normal modes in 


the r-zplane are given by 


a ln 


(=)? + (5). (2.29) 


Cc 
=—— Jk? +k? =— 
fee r Z 2H a : 


0G 


Table 2-1 lists a few values of a_,. 


rae Te 
sans is 


Table 2-1. A few values of a, . 











In our annulus a@=0.10m, so it is seen that Onn is approximately 60 for m= 


1. Similarly, Z, is 0.05 m, so oi is equal to 62.8, or nearly equal to WK 


4 


Z 


Therefore the frequencies of the first cross mode in either the z or r direction are 
approximately equal and have value in the vicinity of 3.3 -3.6 kHz. Our primary 


C 
interest is in frequencies corresponding to L, =A,or f = 7, = 440 Hz. 
eff 


1. Uniform Cross Section Annular Resonator 
a. Undriven case 
As explained above, for low frequencies the annular resonator can be 


modeled as a straight pipe of length L,, with periodic boundary conditions at the 
two end as shown in Fig. 2-2. Assume that there are two traveling waves in pipe. 
One is propagating in the positive x-direction and represented by 

Peace. (2.30) 
The other propagates in the negative x- direction and given by 


P= Be torre) (2.31) 


The total pressure in the pipe is a sum of these two pressures and given by 
P =P, +P. = Ae") + Belo) (2.32) 


These subscripts get used later in a different way use p and n for positive and negative 


instead of 1 and 2. 


L eg ae? TR ory 


<n ann ee 


Ds 


x=0 x=L 


Mp 
He 





Figure 2-2. Annular resonator and equivalent straight resonator. 


Also, the volume velocity U can be obtained from knowledge of the pressure and the 








pipe’s acoustic impedance, Z = = , where S is cross section area of pipe, 
U,=U,-U_= - re airs ae (2.33) 
Periodic boundary conditions are 
Ee veal Oey (2.34) 
and =U 235) 


where L = L,,, . Substituting Eq. (2.32) into (2.34) gives 





em —] 
a Bon (2.36) 
Substituting this expression into Eq. (2.33) and applying Eq. (2.35) yields 
coskL = 1. (2.37) 
This implies that kL = 2m7 Cnt —t0., 1 5 2c... y (2.38) 


or (= = | (2.39) 


It should be pointed out that L,, has slightly different values for different values of 


m. 


b. Driven case 
In addition to knowing the resonance frequencies of the annular 
resonator, we are also interested in the mode shapes, i.e. the acoustic pressure and 
velocity as function of angular position. There is no unique orientation of the modes 
in an undriven, uniform cross section annulus. Either the system must be driven, or 
a nonuniformity in cross section must be introduced to fix the mode shapes. In the 
undriven uniform cross section case, any superposition of two degenerate 
orthogonal modes ( one proportional to cosm@ , the other to sinm@ ) satisfies the 
boundary conditions. In order to find the mode shape, it is necessary to install a 
driver in resonator. Assume that a branch driver is located at x = 0 as shown on Fig 
2-3. The volume velocity injected by the driver is given by 
U 5 = cue (2.40) 





Figure 2-3. Driver installed straight resonator. 


The boundary condition for pressure is given by 
P_,=P.., (2.41) 


x=0 


which is same as that for the undriven system. This conditions leads to 


A(1—e) + BUl-e™) =0. (2.42) 
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The boundary condition for volume velocity now includes an additional term and is given 
by 








Ug, = Opt ,, (2.43) 
which yields the equation 
A(1-e) -B(1+e%) = 2 (2.44) 
Equations. (2.42) and (2.44) can be expressed in matrix form as 
1-e 1-e/* |i A 0 
é : = Pol : 2.45 
1—- e7 FF — 74 ei in S U p ( ) 


This equation has the form MC =U. The value of the coefficients A,B can be 


determined from C=M™'U. Once A and B are known, the pressure and volume 
velocity are found from Eqs. (2.32) and (2.33). 


2: Constricted Annular Resonator. 
a. Frequency vs area ratio (no driver) 
The analysis of the constricted resonator is more complicated due to its 


non uniform cross sectional area. The geometry is shown in Fig. 2-4. 





x=0 x= ] x=L 


Figure 2-4. Constricted resonator with no driver. 


Calculation of the resonance frequency of a constricted resonator is similar to that 


of the uniform cross section resonator. In order to solve this case, the basic equations 


1] 


of pressure and volume velocity in each section have same form as the uniform cross 


section resonator . 
P, = Actin aber = tae (2.46) 
BP, = Cee Die ae (2.47) 


At the junction of two sections, pressure and volume velocity are continuous which gives 


the following equations 


Pea ]) =e =e (2.48) 
Pi(e0)= PG) (2.49) 
Oo Ore (2.50) 
U (x= 0) = 0 (ee (2.51) 


These equations can be combined into a 4 x 4 matrix in the following form. 








pale e JF! _ e — Jk! _ e Ja! 
e 7 aki = e Jal RP e ~Jkal Ded e Jal = On 252 
Sake Stk, G Go4 
l ad 7 Sik, e~JkL Sik, e JkaL D 
| Sk, Sk, 


This matrix has the form of MC = 0. The resonance frequencies are those 


frequencies that satisfy the equation det[M] = 0. 


Unlike the uniform cross section case, the introduction of a constriction 
splits the degeneracy of the modes. The two modes, called the high frequency and 
low frequency modes, correspond to different mode shapes. (In what follows, the 
word “frequency” be will frequency dropped and the modes referred to, simply, as 
“high” and “low”.) The high mode has a pressure antinode centered on the constriction 
whereas the low mode has a velocity antinode centered on the constriction. The 
difference in frequency between the high and low mode increases with increasing 


length or decreasing cross sectional area of the constriction. 
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The dependence of the frequency splitting is shown in Figs. 2.5 a-c. 
The plots show resonance frequency as a function of area ratio for three different 
constriction lengths. The resonator contains air at 20°C and atmosphere pressure. Area 
ratio 1s defined as the ratio of the cross sectional (open) area in the constriction to 
that of the resonator. In the limits of zero area ratio, i.e. when the constriction is 


completely blocked, the high mode frequency corresponds to 4 = (Lig — Leonstriction ) 


and the low mode frequency corresponds to A=2(L,, -L 


constriction 


). These limits are 


appropriate for the two lowest modes of a rigid, rigid resonator of length Ly — Leoncriction « 
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Figure 2-5a. Resonance frequencies with a long constriction (=e 


is 


calculated resonance freq vs area ratio 
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a 
Figure 2-5b. Resonance frequencies with a medium length constriction Ea 


calculated resonance freq vs area ratio 
480 
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L, 
Figure 2-5c. Resonance frequencies with a short constriction 2] 


b. Mode shapes of a constricted resonator 
Calculation of the pressure distribution in a driven, constricted resonator 
is also similar to that for the driven uniform cross section case. For the driven, constricted 


resonator, the resonator is divided into three sections, as shown in Fig. 2-6. The driver 


is located at x = d. 


ee D 


Ba a 


yy Ye, 





se 
ro 

x=0 x= [ 7 x=L 
x=d 


Figure 2-6. Constricted resonator with driver. 


This system gives rise to the following set of equations for pressure and volume velocity in 


the three regions: 





Pa Acer te a ee) (2.53) 

eG Delia aaa (2.54) 

Ea = Ee io *) + Fel oth) | (2.55) 

Cn — i Ae J(@t-k,x) oo Be Matt hx) (2.56) 

U, = 8s Ce J(a@t—kox) ee De ee (2.57) 

U; = > Ee iors) _ —~ ie = Ke haan (2.58) 

and Up =uUpe™. (2.59) 


The three different sections also give rise to six boundary conditions: 
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P(x =0)=P,(x=Z), (2.60) 


U(x =O) — Use (2.61) 
P= 1) = Pee (2.62) 
U (x =D ao are (2.63) 
P,(x =d) =P,(x=d), (2.64) 
and U,(x=d)+y, =U,(x=4). (2.65) 


Application of these boundary conditions to Eqs. (2.53) - (2.59) yields a set of equations 


which can be combined into a matrix equation of the form MC =U. 





! l 0 0 ~ eal _ kal A 6 
] —] 0 0 ~Saky —JkaL Y k, Jka 0 
S; k 1 S ] k ] B 
eu eth! = eu? —e Jel 0 0 C 0 a 
Ge eg az oie! oa edt! 0 0 p| | °° (2.66) 
Sk, Sk, A 0 
0 0 ed ethd _ gated _ oiked ~p,c 
0 0 ead re ae — gated olked i a Up 


Solving the inverse form of Eq. (2.66) gives the value of the coefficients A through F. 
Substitution of these coefficients into Eq. (2.53)-(2.55) gives the mode shapes. 


(1) Mode Selection and Driver Location. Because of 
dissipation, the quality factor (Q) of a particular mode is finite. This means that, in 
general, when the resonator is driven at a frequency in the vicinity of the resonance 
of either the high or low mode, both modes will be excited. In order to measure 
the mode shape of a single mode, it is necessary to excite only that mode as nearly as 
possible. This can be accomplished to an acceptable degree by proper location of 
the driver. If the driver is positioned at a pressure antinode of a particular mode 
and driven at the resonance frequency of that mode, then that mode is the 
predominant one excited. 

Mode selection is illustrated the next sequence of figures which show the 


calculated frequency response and mode shapes for the resonator with a constriction 
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Lag 


length 
ength of 39 


and area ratio 0.7. For Fig. 2-7, the driver is located 180° from the 


center of the constriction which should correspond to a pressure antinode for the 
high mode. Figure 2-7a shows the frequency response as indicated by computed 
magnitude of acoustic pressure measured 45° from constriction. This figure shows that 


there is only a single peak in the frequency response. The peak frequency, 439.0 Hz, 


L 
corresponds to the resonance frequency of the high mode for a ec long constriction 


with a 0.7 area ratio, as shown in Fig. 2-5b. Figure 2-7b shows the calculated mode 
shape when the resonator is driven at 439.0 Hz. It is seen that there is indeed a pressure 
anitnode at both the constriction and driver location. It is also evident that the standing 


wave ratio is very high, another indication of single-mode excitation. 
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Figure 2-7a. Frequency response of the constricted resonator when driver is located 180° 


from the constriction and driven at the high mode resonance frequency for the 


L, 
oy length, 0.7 area ratio constriction. 
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Figure 2-7b. Mode shape when the resonator is driven at the peak frequency as 
determined from Fig. 2-7a. 


Figure 2-8a shows the frequency response when the driver is moved 
to a location 90° from the constriction. This driver location should be near a 
pressure node for the high mode and pressure antinode for the low mode. The 
frequency response again shows only one peak. Comparison with Fig. 2-5b indicates 
that this mode is the low mode. Figure 2-8b shows the mode shape. Itis seen that 
now there is a pressure node at the constriction. Again, the standing wave ratio is 


high indicating single-mode excitation. 
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Figure 2-8a. Frequency response of the constricted resonator when driver is located 90° 


from the constriction and driven at the low mode resonance frequency for the = 


length, 0.7 area ratio constriction. 
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Figure 2-8b. Mode shape when the resonator is driven at the peak frequency as 


determined from Fig. 2-8a. 
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In the cases of the driver located at 180 degrees or 90 degrees from 
the constriction, the driver excites only one mode, either high or low. If the driver is 
placed at a position other than 180 degree or 90 degree, two modes will be excited 


simultaneously. Figure 2-9 shows the frequency response when the driver is located 45 


ee L f a 
degrees from the constriction. We have used the af 0.1 area ratio constriction to 


make the amount of frequency splitting larger. 
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Figure 2-9a. Frequency response when driver is located 45° from the baf 0.1 area 


ratio constriction. 


We can see from Fig. 2-9a that both the high and low modes are excited. In this case, 

they are well separated, in other words, there is no significant overlap. When there 

is little overlap, it is easy to pick up the resonance frequency of each mode. A constriction 
If 

of small area ratio (0.1) and long size ( of 22.5 degree), as in Fig. 2-9a, separates the 


two modes well enough for there to be little overlap. Figure 2.95 shows the mode 


shape when the resonator is driver at the lower peak frequency. This figure shows 
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that even when the driver is at a location which would allow both modes to be 


excited, when the modes are well separated, single mode excitation is still possible. 
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Figure 2-9b. Mode shape when the resonator is driven at the lower of the two 


peak frequencies from Fig. 2-9a. 


ie 
Figure 2-10a shows the frequency response for the tf, 


length, 0.7 area ratio constriction used in earlier figures. In contrast to Fig. 2-9a, the 
modes show significant overlap. This overlap causes two problems. First, it means 
that a driver located at 45° will excite both modes to a significant degree, even 
when driver at one of the peak frequencies. Second, the peak frequencies no longer 
correspond to the resonance frequencies. Overlapping is most significant for large area 
ratio, short constrictions. 

Figure 2-10b shows the mode shape when the resonator is driven 
at the lower peak frequency in Fig. 2-10a. It is seen that the standing ratio is low in 
this case, indicating excitation of more than one mode. The next section addresses a 


technique to determine resonance frequencies in the case of significant overlap. 
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Figure 2-10a. Overlapped wave form of frequency response with driver 45 degrees from a 


ib 
constriction which has a 0.7 area ratio and af, , (5.625 degree) length. 
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Figure 2-10b. Mode shape when the resonator is excited at the lower peak frequency 
in Fig 2-10a. 
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D. POLE - ZERO ANALYSIS 


For understanding the acoustics of an annular resonator, it is important to find 
accurate values of the resonance frequencies. These resonance frequencies are 
independent of the driver location but are dependent upon what kind of constriction is 
being used - the length and area ratio. As we have seen, the driver may excite one or 
two modes simultaneously depending upon its excitation frequency and location. When 
only one mode is excited, it is simple to find resonance frequency by finding the 
frequency corresponding to the peak in the frequency response . But for the case when 
the driver excites two modes simultaneously, the peak frequency cannot be used to 
determine the resonance frequency. This situation was discussed in the previous section. 
Pole-zero analysis can be used to determine the resonance frequencies in the case of 


significant modal overlap. 
i Analysis 


For a driven simple harmonic oscillator, the complex speed is represented by [Ref 
1] 


jat 
ae eee ee (2.67) 


R + if s) 
ni--—_ 
os J | @ > 


The expression for 7 can be written as 


: He jot 
- ja —e 
= m 
“u= Sy a (2.68) 
a - jo -@O; 
Q 


by introducing the resonance frequency w, and quality factor, O = oo 


By analogy, we assume the pressure in a resonator driven near resonance, at 


One point in space, has the form 


Zs 


2 —————— (2.69) 


@ 
o* -~j—-o-o: 


oO 


Q 


where Q is now the quality factor for the resonator. Factoring the denominator of Eq. 


(2.69) results in 





A 
2 =. (2.70) 
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This equation can be simplified if the denominator is written as [Ref 2] 














(@ + ip, Xo + jP.), (2:71) 
where the poles p, and p, are, 
2 
py “met j.|@, (2) (2.72) 
2 
and P= + j,{o; - (2 (273) 
2 
or p=- sat j,lo? - (2 (2.74) 


where @, = @,, for the single mode case. The poles are written in the form (@ + jp ) 


because the Hewlett Packard signal analyzer used for experimental measurement outputs 


the poles in that form. Equation (2.74) can be simplified into a form given by 


2 
p=-at jb where @ = aan and b = ao, = (o ; (2.75) 
] ] 
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For two overlapping modes, the complex pressure amplitude is the sum of the individual 


amplitudes 
= P,+P, (2.76) 
and P,, has the form 


_ Ai(@ + JP3)(@ + Jpg) + Az (@ + JP: @ + JP2) 


P 
(@ + Jp, > (@ + jp,) 


(Qaim 


tot 


So, even though the two modes are overlapped and the amplitude of the individual 


responses are added, the poles still show separable frequencies. 


The numerator of Eq. (2.77) can be factored in the same way as the denominator 


of Eq. (2.69). In this case the numerator is 


eset lee} 


79 
0. ), (2.79) 


o 
where c= A (-—)+A,(- 
2 








and d= -)a(=$4) +4,(34)| +4(A,+A,\A,o2+A,07). (2.80) 


From Eq. (2.78) we can find zeros. 


The resonance frequency and quality factor can be obtained from the poles using 


Eq. (2.75). The resonance frequency is then 


QO, = ala? = b? (2.81) 


2» 


and the quality factor is 


Q,=—. (2.82) 


Q, 


20, 





If Q is large, an approximation can be used to simplify Eq. (2.81). The term of 





2 

‘ @ : 

can be neglected since ( 5 ' is small compared tow”. So, b = @,. 
1 

For data analysis, we fit the measured frequency response using pole-zero 


analysis. The fit parameters give the resonance frequency and Q of the individual modes. 
E. THERMO-VISCOUS DAMPING 


Up to this point, losses have been ignored. However any real gas has viscosity 
and will exchange heat with the walls of the resonator. The boundary condition at 
the walls of the resonator is that the normal velocity is zero. When viscosity is included 
the correct boundary condition has both the normal and tangential velocity going 
to zero. This creates a viscous boundary layer which contributes viscous losses to the 
propagation of sound [{Ref. 3,4]. Similarly, the gas near the walls will exchange 
heat with the walls. This thermal boundary layer also contributes losses to the sound 
propagation [Ref. 3,4]. While these viscous and thermal effects are complicated, if the 
boundary layer thicknesses are small compared to the cross sectional dimensions of the 
resonator, the effects can be included by a modification of the wave number k. 


The frequency dependent viscous and thermal boundary layer thickness 6, and 


O, are given by 


2u 
= |-—— 2.83 
6, =e (2.83) 
Pp 
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where yw is the viscosity, K is the thermal conductivity, c,is the specific heat and p 


is the density of the gas. The wave number k is then given by [Ref. 5] 


a2 2m +7 D4 J) 
: c 2r, ~6,(1+ J) cS 


where r, is the hydraulic radius, which is the ratio of the cross sectional area to the 
perimeter. It should be noted that since the hydraulic radius is different in the constriction 


and in the main channel, and that k is now complex. The imaginary part of k is related 


to the thermo-viscous damping. 
F. END EFFECTS 


In the earlier analysis, the boundary conditions used at the ends of the constriction 
were continuity of pressure and volume velocity. Using only plane wave analysis, the 
pressure and particle velocity are constant across the cross section. However, at the 
constriction end, the velocity and pressure have a spatial dependence and thus the plane 
wave theory is not as accurate near the ends as it is away from them. 

The spatial variation in the pressure can be described using non-propagating 
higher order modes. These non-propagating modes store energy near the constriction 
ends. Their effect on the plane wave propagation can be modeled as a lumped impedance 
element. [Ref. 6]. The lumped impedance can be obtained from a variety of methods. 
A first approximation involves modeling the impedance as an added mass much like 
the end correction for radiation at the end of a tube. This added mass can be computed 
using conformal mapping. [Ref. 6]. A more accurate method involves including higher 
order modes near the constriction end. [Ref. 7]. While computationally more intensive 
than the conformal mapping solution, the higher order mode solution is more accurate. 

After computing the lumped impedance Z(@), the boundary condition for the 
pressure becomes (the volume velocity boundary condition is unchanged.) 

P, (J) = P,() + Z(@)U,(/) (2.86) 


as 


where P, is the pressure in unconstricted region and P, is the pressure in constriction. 


U,, is the velocity in the constriction. The constriction is at x = /. 
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Hil EXPERIMENTAL APPARATUS AND PROCEDURE 


A. INTRODUCTION 


A description of the experimental apparatus and the procedures for data 
acquisition are presented in this chapter. The first section will include a description of the 
annular resonator, constrictions, the drver, the microphone and_ electronic 
instrumentation. The second section discusses the data acquisition system. The third 


section is a detailed discussion of the experimental procedure. 


B. EXPERIMENTAL APPARATUS 


1. Resonator and Driver 


Figure 3-1 is an illustration of the annular resonator and acoustic driver. The sides 
and bottom of resonator are made of aluminum. The top of resonator is made of 
plexiglas to allow visual sighting of the relative positioning of the constriction and the 


driver. More details about each part will follow. 


The dimensions of the resonator were chosen to give longitudinal resonance 
frequencies in the 300 ~ 500 Hz range and to be well below the cutoff frequency 
for cross modes. The resonator consists of two main parts as shown on Fig. 3-1. 
One piece is an aluminum ring with 15.30+0.02 cm inner radius and 1cm thick wall. 
This ring formed the outer side wall of the resonator. The other piece is a solid 
aluminum cylinder which is concentric with the outer ring. The inner radius of 
resonator is 10.040.01 cm. The height of resonator is 5.040.01 cm. The resonator 
is capped on the open side with lcm thick plexiglas plate with a radius of 16.0 cm. 
The bottom plate of the resonator is made of aluminum and its thickness is 1.0 cm. 
Sixteen microphones are installed in the cap to measure the spatial acoustic pressure 


distribution inside resonator. The installation of microphones will be discussed later. 


Zo 


The cap is bolted with 10.7cm long threaded rod in the center and 8 clamp bars as 
shown on Fig. 3-2. To ensure a tight seal (for a high quality factor), vacuum grease 


was applied to the surface of the aluminum ring and center piece. 


View of resonator from above without cap 


Threaded rod 


TTT 


is 


Extra clamp bar 





Figure 3-1. Llustration of resonator and driver. 


The acoustic signal was generated by a 40 W , 16 ohm external driver. The 
driver was connected to the resonator by a 3.5 cm long piece of plastic hose with 0.4 
cm diameter. The driver was excited by the output of a Hewlett Packard 3562A dynamic 
signal analyzer sent through a Techron 5507 power amplifier. The signal analyzer 
was used to determine the frequency response of the resonator. Typically, the 
resonator was driven with random noise with a 1 Vrms source level. To measure the 
mode shape, a Hewlett Packard 33120A function generator was used to excite the 
resonator at the resonance frequency determined from the built-in pole-zero analysis 


of the frequency response using the signal analyzer. 


30 


Clamp bar Threaded Rod & nut 
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Figure 3-2. Illustration of clamp bar, threaded rod and nut. 





2. Microphones 


For the experiment, 16 Panasonic WM60-AY microphones were used to 
measure the pressure distribution inside resonator. The dimensions of the microphones 
are 6mm in height and 5mm in diameter. The sensitivity is -42 dB + 3dBre1Pa. A 
Hewlett Packard 6237B triple output power supplier provided 9 V DC to each 


microphone. 


a. Installation of microphones 

The microphones were installed as shown on Fig. 3-3. The first 
intention was to use only one microphone and rotate the cap for each different position 
to measure the mode shape. However, the rotation of cap broke the sealing between the 
plexiglas cap and aluminum bodies of the resonator making the measurements very time 


consuming and nonrepeatable. 
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Figure 3-3. Installation of 16 microphones on the plexiglas cap. 


We decided to install 16 microphones in order to protect sealing and reduce the time 
required for data acquisition. For the precise measurement in constricted area, 6 
microphones are located in first 30 degree region where a constriction is located. The 
remaining 10 microphones are located 30 degrees apart. The microphones were installed 
in the center of the annular channel. To install the microphones, 16 holes were drilled 
part way through the plexiglas cap. These holes were slightly larger than the diameter of 
the microphone. A smaller hole was drilled concentric with each of the larger holes. This 
smaller hole penetrated the cap. Therefore the hole size presented to the inside of the 
resonator was much smaller than the microphone diameter. The small hole diameter is 
0.15 cm. The microphones were sealed into the holes with silicon glue. The error in 


the microphone position is estimated to be less than + 0.25°. 
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Figure 3-4 . Chamber used for microphone calibration 


b. Calibration of microphones 

Because 16 microphones are used for sensing the pressure distribution, the 
first step was to calibrate all the microphones over the 300 ~ 500 Hz frequency range 
interest. For the calibration, a custom made chamber with a built in speaker was 
used. As shown on Fig. 3-4, the diameter of the chamber is 3.0 cm and the height is 2.50 
cm. At first we tried to measure the pressure for each microphone at the same 
position relative to the driver in an otherwise empty resonator. However, the phase 
shift caused by the resonators’ resonance made calibration very sensitive to small 
errors in position. Instead, the top cap has removed and the calibration chamber 
was positioned directly over each microphone. The frequency response was measured 
for each microphone over a 300 ~ 500 Hz range. A typical frequency response is 
shown in Fig. 3-5. This sweeping and data storage were done by Hewlett Packard 
35665A dynamic signal analyzer. For the calibration, the signal analyzer was set to swept 
sine mode with a settling time of 20 cycles and integration time of 50 cycles. The 
chamber was placed onto the cap where each microphone is located and the chamber 


was sealed on the cap with vacuum grease. 
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© X:353 Hz Y:474.0161 m 





Figure 3-5. Frequency response of Mic. #1 covered by the calibration chamber. 


Microphone #1 was used as the reference for comparison calibrations of the 
other 15 microphones. Each microphone’s output was divided by that of 
microphone #1 to create a correction factor as shown on Fig. 3-6. When measuring 
the mode shape, the output from each microphone was divided by the appropriate 
correction factor to obtain the corrected amplitude at a particular location and 
driving frequency. For example, consider the case when the measured amplitude from 
microphone #12 is 50 mV at 455 Hz. The correction factor is 0.739 from the Fig. 
3-6. Dividing 50 mV by 0.739 gives 67.66 mV as the corrected measurement. 


Generally, the correction factors were between 0.7 and 1.2. 


© X:455 Hz Y:738.7183 m 
FRES/D1 2:1 





Figure 3-6. Plot of correction factor for Mic. #12 referenced to Mic. #1. The correction 


factor is 0.739 at 455 Hz. 
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3: Constriction 


To explore a wide range of parameters, various sizes of constriction are necessary. 
When considering the material selection for this constriction, a couple of factors came 
to mind. The material had to be rigid and have good machining propertied for ease 
of fabrication. These two factors led to the selection of PVC (Polyvinylchlonde). 
Figure 3-7 shows the design along with the nominal and measured dimensions for the 
various constrictions used. For the experiment, 3 different angular widths were used. 


The width of constriction is given by the @ shown in Fig. 3-7. The three widths 


are 22.5° , 11.25° and 5.625°. The 22.5° constriction corresponds to — of the 
effective circumference of the resonator. 11.25° corresponds to = of L,, and 


5.625° corresponds to = of L,,. The area ratio is the ratio of open area in the 


constriction (which is the shaded black square in Fig. 3-7) to cross sectional area of 
resonator. In Fig 3-7, $1 stands for the open area and S2 stands for the open cross 
sectional area of the constriction. The table in Fig. 3-7 gives the geometrical parameters 


of 32 different constrictions. 
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unit: cm 
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Figure 3-7b. Measured geometry of the constrictions. 


4, Electronic Instrumentation 


Figure 3-8 is a block diagram of the electronic instrumentation used during 
the measurements. A Hewlett Packard 3562A Dynamic signal analyzer is used to 
determine the frequency response in the band from 300 Hz to 500 Hz. The signal 
from this analyzer is amplified by a Techron 5507 power supply amplifier and used as 
the input signal to the driver. The pole-zero analysis function built into the signal 
analyzer is used to extract the resonance frequencies from the measured frequency 
response. Once the resonance frequency is found, a Hewlett Packard 33120A 
Function/Arbitrary Waveform Generator is used to generate a 1 Vpp_ sinusoidal 
signal at the obtained resonance frequency. The signal is also amplified by the 
same Techron amplifier and used as the input signal to the driver. Either the signal 


from the analyzer or the function generator signal is monitored on a Tektronics 2336 


on 


Oscilloscope. The outputs from the microphones are monitored on the same oscilloscope. 
A custom made microphone switch board selects each microphone to measure the 
pressure amplitude at different positions. The microphone output amplitude is measured 
using a Standford Research 530 Lock In amplifier. In order to measure the 
temperature of the resonator, an Omega 450 AET thermocouple thermometer type E 1s 
used. The thermocouple was calibrated to +0.1°C using a distilled water ice bath. The 
thermocouple is positioned on the surface of the outer aluminum ring of resonator. 
The temperature inside resonator is measured frequently and be recorded whenever there 


iS variation. 
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c EXPERIMENTAL PROCEDURE 
rr Resonator Setup 


Prior to any data acquisition, the resonator is assembled with the constriction in 
the desired position relative to the driver. For the experiment, three major positions are 
used. They are 45°, 90° and 180° from driver. The position is measured from the 
driver to the center of constriction. The top cap is oriented so that 0° is set at the 
left edge of constriction. The 0° setting is performed by rotation of the cap by hand. 
For the error measurement of 0° setting, three people conducted rotation of the 
cap. Each person tried three times and we found that the error by rotation was less 
than +0.5°. So, with the mechanical error +0.25° in dniling the microphone hole, the 
total error of microphone positioning is less than +0.75°. With the constriction 
positioned, the proper amount of vacuum grease is put on the open surface of resonator 


and the top capis secured with the clamp bars. 
2. Determining Resonator Characteristics 


The Hewlett Packard 3562A Dynamic Signal Analyzer is used to characterize 
the resonator. The swept sine mode of the analyzer provides a sine wave signal to 
the system at a given range of frequencies. The frequency response of the system is 
used to determine the resonance frequency and quality factor of several modes of the 
resonator. When only one mode is excited, the peak pressure amplitude frequency is 
same as the resonance frequency given by pole-zero analysis. When two modes are 
excited, pole-zero analysis is used to determine the resonance frequency. The signal 


analyzer gives polesin the form a+ jb. The resonance frequency is simply the 


b : 
value of 5. The quality factor can be also obtained by Bn from that form. This 


method has been explained in the theory chapter. These measurements are performed 


for each constriction type and each relative position of constriction and driver. 
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3. Data Collection 


Data collection is performed insets of the positions and different sizes of 
constriction. Before data collection, the initial ambient temperature and noise level 
is checked. The output of the Hewlett Packard 6237B Triple output power 
supply is also set to DC 9V. The measurement is then commenced by getting the 
resonance frequency from the Hewlett Packard 3562A Dynamic signal analyzer. The 
microphone switch board is set to the microphone which is closest to the driver. Then 
the resonance frequency and the temperature of resonator was measured. The 
procedure to acquire resonance frequency has been mentioned in subsection 2 in 
this section. Getting quality factor also follows the same procedure previously 
mentioned. After finding resonance frequency, a data run is commenced by 
hooking up the Hewlett Packard 33120A function generator to the input amplifier for 
driver and setting the frequency and voltage. The driving frequency should be the 
resonance frequency obtained by previous procedure and the voltage is 1 Vpp. The 
acoustic amplitude is measured with each microphone. With the driver still running, the 
switch board selects a particular microphone and the value of amplitude at each 
position is measured with the SR 530 Lock - In amplifier. At this time, if there is 
variation in the resonator temperature, it is noted. All data recorded from the 16 
microphones are scaled by the correction factor obtained from microphone calibration. 
The corrected data go to computer program to plot mode shape. This concludes data 


acquisition of one mode. 
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Procedure for Wave Guide Characterization Collection 


1. Date: ; [97 ( ~ ) 

2. Ambient Condition 
Room Temper : [Ct K) / Speed of Sound = m/s 
Ambient Pressure: mbar / Initial Temp of resonator: Ae 


3. Configuration (Constriction) 


# of Const. : EA / Size : e 
Area Ratio : 
Position of Const. : ° (CCW/CW) from driver 
Position of Driver : 5 
4. Resonance Freq: / Q: / Mode: 
5. Source information 
Input : Vpp / Type of source : 
6. Data: WG with Const(s mode 
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Figure 3-9. Data collection form. 
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IV. RESULTS AND DISCUSSION 


The purpose of this chapter is to show the comparison between theory and 
experimental results. The result of resonance frequency measurements are discussed first 


and then the mode shapes. 


A. RESONANCE FREQUENCIES 


The measured resonance frequencies are determined from pole-zero analysis of the 
frequency response. Because the resonance frequencies were measured at different 
temperatures, they were corrected to 20°C in order to make comparison of data 
recorded at various temperatures with the theoretical resonance frequencies more 


straightforward. The equation for conversion is 


ro eee | aa (4.1) 





The theoretical resonance frequencies were calculated for two different cases. One case 
was without end corrections and the other was with end corrections. All the results are 
listed in Tables 4-1 and 4-2. Tables 4-la and 4.2a list the resonance frequencies for the 
low (L) and high (H) modes for the various constriction lengths and area ratios. Tables 4- 
1b and 4-2b lists the temperatures at which the frequencies were measured. Tables 4-1c 
and 4-2c list the 20°C-corrected frequencies. Tables a, b and c are identical for Tables 4-1 
and 4-2. Tables 4-1d and 4-2d give the calculated resonance frequencies. Table 4-1d 
ignores end corrections. Table 4-2d includes end corrections as calculated from higher 
order non propagating modes. Tables 4-le and 4-2e list the percent error between 
measured and calculated resonance frequencies. 

It is evident that the inclusion of end corrections has little effect on the prediction 
for the high modes. In both cases the average error is approximately 0.15%. However, 
end corrections are very important for the low mode. This point is evident from Fig. 4-1 


where it is seen that inclusion of end corrections is important for the low mode only. 


43 


Figure. 4-1 shows the 20°C-corrected resonance frequencies as well the predicted results 
(without and with end corrections) as functions of area ratio. 

In this section it was shown that end corrections are an important factor to include 
when calculating resonance frequencies, especially for the low mode. In the next section, 


it is shown that they are also important when calculating mode shapes. 


i Comparison of Measured Resonance Frequency with Calculated Results 


without End Corrections 






TART 01] 02] 03] 05 | 07 10 
eae eH | H | Liao Lt i an | 
| 1/32 _| 328.7 | 449.6 | 368.7 | 447.8 | 391.7 | 446.4 | 416.5 | 443.6 | 426.9 | 439.3| 
1/64 | 348.6 | 443.1 | 384.2 | 441.4 | 402.6 | 440.7 | 422.1 | 439.3 431.1[ 438.6| 


Table 4-1a. Measured resonance frequencies for different constrictions. 











AR [on | 02 oS 0 
size | LHL] a] te} AT tH] Lh | a 
Pe ee (227s 23/3 23652316 23-9 23.7) 21.60) 21-6) | 
240 2400922322 301203 22.3223 |) 22.3) 23:59] 23.5 


Table 4-1b. Temperatures for each measurement from table 4-1a. 












a 
L/32 | 326.5 | 446.6 | 366.6 | 445.3 | 389.3 | 443.7 | 413.9 | 440.8 | 425.7[ 438.1 | 
L/64 | 346.2 { 440.1 | 382.7 | 439.7 | 401.0 | 439.0 | 420.5 | 437.6 | 428.5| 436.0] 


a = Measured resonance frequencies corrected to 20°C. 











ee LA H | _L 
1/32 | 349.3 | 445.9 | 387.3 | 444.5 | 404.2 | 443.1 | 420.2 | 440.2 [427.1 [4377| 
1/64 [381.1 | 439.5 | 407.7 | 438.8 | 417.8 | 438.1 | 426.5 | 436.7 | 430.1] 435.5) 


Table 41d. Calculated resonance a ee at 20°C. ato aa corrections. 
















PART oor To? ts 
psie | L {| Ht ot | wf fe 
| ue | 3.48 [0.07 | 3.71 [0.18 | 2.91 | 0.21 | 1.29 | 0.13 | 0.32 | 0.04 | 0.05 
(132 [698 | 0.16 | 5.65 | 0.19 | 383 | 0.14 [153] 013 | 033010; 
[1/64 | 10.1 | 0.13 | 6.54 | 0.20 | 4.19 | 0.20 [1.44 | 0.20 [038 [ou] 


- 4-le. Percent error between measured and calculated resonance frequency. 
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ze Comparison of Measured Resonance Frequency with Calculated Results with 


End Corrections 


| 1/32 _| 328.7 | 449.6 | 368.7 | 447.8 | 391.7 | 446.4 | 416.5 | 443.6 | 426.9 | 439.3 | 
_1/64_| 348.6 | 443.1 | 384.2 | 441.4 | 402.6 | 440.7 | 422.1 | 439.3 | 431.1] 438.6 | 


Table 4-2a. Measured resonance frequencies for different constrictions. 















AR | 01 | 02 

Tse { L[A|L[A|L[]A|[LE[AI]L [Al 
1/32 | 23,97) 2319)| 23:35) 23:3) 12316) 236 )23:7 ase ee a 
_ 1/64 | 24.0 || 24.0 | 22:37} 22/3) |) 22131)))2213))))22:30)223) 3502S 


Table 4-2b. Temperatures for each measurement from table 4-2a. 












PAR | 01) | 02 | 0ST 
ee ee 
1/32 | 326.5 | 446.6 | 366.6 | 445.3 | 389.3 | 443.7 | 413.9 | 440.8 | 425.7] 438.1 | 
_L/64 | 346.2 | 440.1 | 382.7 | 439.7 | 401.0 | 439.0 | 420.5 | 437.6 | 428.5] 436.0] 


Table 4-2c. Measured resonance frequencies corrected to 20°C. 


ris2 [3267 | 445.9] 3662 | asda | 3887] 443.0] 4141] 4402 [209 [@a7.7[ 
| 1/64 | 345.6 | 439.5 | 381.4 | 438.8 | 400.0 | 438.1 | 420.0 | 436.8 | 427.81 435.5] 


Table 4-2d. Calculated resonance frequencies at 20°C with end corrections. 













_sie | L | A |b | ao) asa eee ae 
Lis | 034 | 0.08 | 028 | 021) 020) 02) 
(132 [0.10 [0.16 | 0.11 | 0.19 | 0.15 | 0.15 | 0.04 | 0.14 | 018 | 0.09 | 
_Le4 | 0.17 | 0.13 | 0.34 | 0.20 | 0.24 [0.20 [ 0.12 | 019 [015 | 012 [ 


Table 4-2e. Percent error between measured and calculated resonance frequency. 
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measured & calculated resonance freq vs area ratio 
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Figure 4-la. Comparison of calculated and measured resonance frequencies versus area 


ratio for the L/16 constriction. 


measured & calculated resonance freq vs area ratio 
480 


460 
440 
420 
400 

f 380 
360 
340 


320 





soa 0.2 0.4 0.6 0.8 1 


area ratio/const:L/32 


Figure 4-1b. Comparison of calculated and measured resonance frequencies versus area 


ratio for the L/32 constriction. 
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measured & calculated resonance freq vs area ratio 


480 
460 
440 
420 
400 

f 380 
360 
340 


320 





300 
0 0.2 0.4 0.6 0.8 1 


area ratio/const:L/64 


Figure 4-lc. Comparison of calculated and measured resonance frequencies versus area 


ratio for the L/64 constriction. 


B. MODE SHAPES 


Figures 4-2a and b show the comparison of the measured and computed mode 
shapes for the L/16, 0.1 area ratio constriction. The measurements have been normalized 
to the highest value of pressure amplitude. The amplitude of the mode shape is 
determined from least squares fit to the measured data. The driver is located 45° from the 
constriction. Figure 4-2a shows the results for the high mode. It is seen that there is 
good agreement between the measured and predicted mode shapes, even though end 
corrections are not included in the calculation. The driver is located at a place where both 
modes should be excited. However, the modes are well separated with this constriction, 
so there is little excitation of the low mode. Single mode excitation is indicated by a high 
standing wave ratio from the mode shape. 

Figure 4-2b shows the mode shape for the low mode. The agreement between the 


measured and calculated mode shape is significantly worse than with the high mode. 
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Figures 4-2c and d show the mode shapes again but this time end corrections are 
taken into account. Comparing Fig. 4-2c and 4-2a, it is seen that end corrections have 
little effect on the high mode. However, they are very important for the low mode, just 
as was discussed in the previous section. Figure 4-2d shows that when end corrections 
are taken into account, the agreement between calculated and measured mode shape 
improves significantly. 

Figures 4-3a and b show the mode shapes for the case when the driver is located 
45° from the L/32, 0.7 area ratio constriction. In the case the constriction is relatively 
short and the area ratio is relatively large. As a result, the modes are not well separated. 
So although the driver is excited at the resonance frequency of one mode, there is 
significant excitation of the other mode, as evidenced by the low standing wave ratio. 

Figure 4-3a shows the mode shape when the resonator is driven at the high mode 
resonance frequency. End correction are ignored in the calculation. The agreement 
between measured and calculated mode shape is not as good as in the previous, well 
separated-mode case (Fig. 4-2). Figure 4-3b shows the results for the low mode. Here 
the agreement between theory and experiment is worse. In both cases, the standing wave 
ratio is noticeably lower that in the well-separated case. 

Figure 4-3c and d show the effect of including end corrections. There is some 
improvement in the agreement for the high mode case. There is considerable improvement 
in the agreement for the low mode, and some improvement for the high mode. Many other 
comparisons are provided in the Appendix. The cases presented in this chapter are 
representative. The overall agreement in mode shape is good even in the overlapping- 
mode case. 

The results of this chapter indicate that end corrections are important both for 
calculation of resonance frequency and mode shapes. Also end corrections are more 
important for the low mode. This is because the acoustic impedance of the low mode is 
low in comparison to that of the high mode. The addition of the effective impedance of 


the end corrections is of more consequence, then, to the low mode. 
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Mode Shape for Length L/16, Area Ratio 0.1, and f=461.5 


— Theory 

*% Measurement 
Edge of Const. 

-+- Driver Position 
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Theta (degrees) 


Figure 4-2a. Mode shape for the high mode when the driver is located 45° from the 
constriction having length of L/16 and area ratio of 0.1. End corrections are not included 


in the calculation. 


Mode Shape for Length L/16, Area Ratio 0.1, and f=308.2 
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Figure 4-2b. Mode shape for the low mode when the driver is located 45° from the 
constriction having length of L/16 and area ratio of 0.1. End corrections are not included 


in the calculation. 
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Mode Shape for Length L/16, Area Ratio 0.1, and f=461.5 
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Figure 4-2c. Mode shape for the high mode when the driver is located 45° from the 
constriction having length of L/16 and area ratio of 0.1. End corrections are included in 


the calculation. 
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Figure 4-2d. Mode shape for the low mode when the driver is located 45° from the 


constriction having length of L/16 and area ratio of 0.1. End corrections are included. 
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Mode Shape for Length L/32, Area Ratio 0.7, and f=439.3 
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Figure 4-3a. Mode shape for the high mode when the driver is located 45° from the 


constriction having length of L/32 and area ratio of 0.7. End corrections are not included. 


Mode Shape for Length L/32, Area Ratio 0.7, and f=426.9 
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Figure 4-3b. Mode shape for the low mode when the driver is located 45° from the 


constriction having length of L/32 and area ratio of 0.7. End corrections are not included. 
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Mode Shape for Length L/32, Area Ratio 0.7, and f=439.3 
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Figure 4-3c. Mode shape for the high mode when the driver is located 45° from the 


constriction having length of L/32 and area ratio of 0.7. End corrections are included 


Mode Shape for Length L/32, Area Ratio 0.7, and f=426.9 
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Figure 4-3d. Mode shape for the low mode when the driver is located 45° from the 


constriction having length of L/32 and area ratio of 0.7. End corrections are included 


53 





ve SUMMARY AND CONCLUSIONS 


The purpose of this thesis is to investigate frequency splitting and mode shapes in a 
constricted annular resonator. The measurements are intended to serve as preliminary 
verification of a model predicting the performance of an annular thermoacoustic prime 
mover. 

The constriction splits the degeneracy of the plane wave-like modes of a uniform 
cross section annular resonator. One mode has a pressure antinode at the constriction, the 
other a pressure node. The resonance frequency of the former is higher than that of the 
latter. The difference in resonance frequency depends upon the length and porosity of the 
constriction. The effects of dissipation cause the frequency responses of the individual 
modes to overlap. 

An equivalent plane wave model is developed to predict the resonance frequencies 
and mode shapes. The experimental resonance frequencies are determined from pole-zero 
analysis of the measured frequency response. 

It 1s shown that it is critical to take into account end corrections at the constriction 
to get good agreement between the measured and predicted results. End corrections are 
more important for the low mode than the high mode. 

In general the percent error between measured and predicted resonance 
frequencies is less than 0.2 % for high mode and 0.3 % for the low mode. The measured 
mode shape are also in good agreement with predictions. 

The results of this thesis indicate that the simple constricted annular resonator 1s 
well understood. The sensitivity of the resonance frequencies and mode shapes to end 
corrections indicates that the annular resonator is a good geometry for investigating these 
effects further. For instance, how does one model the end corrections for a thermoacoustic 
stack? Are the end corrections for a stack and a constriction having the same porosity 


substantially different? These questions are left for future work. 
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APPENDIX DATA TABLE & MODE SHAPE 


Experiment 1. 














, arearatio:0.1 , ang e from driver:180° 


Temp: 24.6°C 
Sele ee ete 
# (6 (mV) Data # (0 (mV) Data 
10 [705 [705 | 9 a2) | 31 | 166 
3 @a) | 43 | 4 | Gm | 31 | 80 
[6 3) | 389 | ol | 14 em | 58 | 77 





Mode Shape for Length L/16, Area Ratio 0.1, and f=463.7 


— Theory 
* Measurement 
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Experiment 2. 















h:L/ 16, arearatio: 0.1 , angle from dnver:90° 


# (6 (mV) Data # (0 (mV) Data 
2 (06) | 172_| 176 | 10 so) | 198 | 295 
7060) | 449 | 606 | 15 G00) | 349 | 578 


Mode Shape for Length L/16, Area Ratio 0.1, and f=308.8 
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Experiment 3. 


length: L/16 , arearatio:0.1 , angle from driver:45° 


Temp: 22.5°C 
# (0 (mV) Data # (0 (mV) Data 


Mode Shape for Length L/16, Area Ratio 0.1, and f=461.5 
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Experiment 4. 


















length: L/ 16, arearatio: 0.1 , angle from driver:45° 


Temp: 22.5°C 
# (6 (mV) Data # (0 (mV) Data 
4 | 77_| 30 | 2 Gi | 89 | 127 


Mode Shape for Length L/16, Area Ratio 0.1, and f=308.2 
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Experiment 5. 


, arearatio: 0.1 , angle from driver: 180° 


# (0 (mV) Data # (0 (mV) Data 
1 00) | 96 | 6 | 9 (120) | 206 | 260 
3) | 352_| 705 | 1 (a | 28 | 099 
805) [76 | 99 | 16 G30) | 502 | 590 


Mode Shape for Length L/32, Area Ratio 0.1, and f=448.9 
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Experiment 6. 


length: L/32 , arearatio: 0.1, angle from driver:90° 


# (0 (mV) Data # (6 (mV) Data 


502 
5 G30) 
7060) 
500) 


Mode Shape for Length L/32, Area Ratio 0.1, and f=328.7 
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Experiment 7. 
















length: L/32 , arearatio: 0.1 , angle from driver:45° 


# (6 (mV) Data # (0 (mV) Data 
[1 oo [447 [a7 [9 2) | 133 [68 
[7 (050) [210 | 284 [15 Goo) | 199 | 206 


Mode Shape for Length L/32, Area Ratio 0.1, and f=449.6 
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Experiment 8. 















length: L/32  , arearatio: 0.1, angle from driver:45° 


Temp: 23.8°C 
# (0 (mV) Data # (6 (mV) Data 


Mode Shape for Length L/32, Area Ratio 0.1, and f=328.7 
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Experiment 9. 


, arearatio: 0.1, angle from driver:180° 


# (0 (mV) Data # (0 (mV) Data 
2 (0s |_o77_[ 099 | 10 (150) | #14 | 617 


& 3 | 69 | 066 | 14 G7) | 14 | 18 
7060) 341 
30%) 2 


Mode Shape for Length L/64, Area Ratio 0.1, and f=443.2 





IPI 


— Theory 

* Measurement 
- -- Edge of Const. 
--« Driver Position 





0 50 100 150 200 250 300 350 
Theta (degrees) 


65 


Experiment 10. 


Temp: 24.0°C 
Se Se Ste 
# (6 (mV) Data # (0 (mV) Data 


, arearatio: 0.1, angle from driver:90° 


630) | 212 [248 | 14 G70) | 198 | 264 


Mode Shape for Length L/64, Area Ratio 0.1, and f=350.2 
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Experiment 11. 


length: L/64_ , arearatio: 0.1, angle from driver:45° 


# (0 (mV) Data # (6 (mV) Data 
sy [39 [a4] 3 toy te | 209 
sa 310] 363.4 rm) [08] 
7060) | 168 | 7 | 15 G00) | 186] 193 


Mode Shape for Length L/64, Area Ratio 0.1, and f=443.1 
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Experiment 12. 


length: L/64 , arearatio: 0.1 , angle from driver:45° 


# (0 (mV) Data # (0 (mV) Data 
2008) [127 | 13.0 | 10 (is) | 128 | 190 


3) [205 [259 | 1 Gia) | 21 | 8 





Mode Shape for Length L/64, Area Ratio 0.1, and f=348.6 
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Experiment 13. 
length: L/ 16 


9 , arearatio: 0.7 , angle from driver:180° 


# (6 (mV) Data # (0 (mV) Data 
2 (wey |_o76 [| 9 | 1050) | 347 | si7 
[4 crs |_331_| 70s | 12 Gio | «6 | 659 
sy] 4s | 6 | Bem | 54] 59 


Mode Shape for Length L/16, Area Ratio 0.7, and f=444.2 
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Experiment 14. 
length: L/ 16 


Q 
Temp: 22.1°C 
# (6 (mV) Data # (0 (mV) Data 
401g) [92 | a3] 12-1) | 166 | 237 
sm] 02139] 3 aay | aa | 5 


Mode Shape for Length L/16, Area Ratio 0.7, and f=421.2 


, arearatio: 0.7 , angle from driver:90° 
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Experiment 15. 


, arearatio: 0.7 , angle from driver:45° 


Se Sle se 
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Mode Shape for Length L/16, Area Ratio 0.7, and f=444.0 
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Experiment 16. 














, arearatio: 0.7 , angle from driver:45° 


Sele ete se 
# (0 (mV) Data # (0 (mV) Data 
[205 | 67 | 69 | 10 50) | 258 | 384 


Mode Shape for Length L/16, Area Ratio 0.7, and f=421.9 
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Experiment 17. 


, area ratio: 0.7 , angle from driver:180° 


Se ee 
# (6 (mV) Data # (0 (mV) Data 
[2 (00 | 74 | 096 | 10 (130) | 375 | 559 


Mode Shape for Length L/32, Area Ratio 0.7, and f=441.1 
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Experiment 18. 














length: L/32 , arearatio:0.7 , ang 


Temp: 21.1°C 
# (0 (mV) Data # (0 (mV) Data 
302) | 107 | 6 | mcs | 72 | 61 


Mode Shape for Length L/32, Area Ratio 0.7, and f=426.2 
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Experiment 19. 
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Mode Shape for Length L/32, Area Ratio 0.7, and f=439.3 
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Experiment 20. 














, arearatio: 0.7 , angle from driver:45° 


Temp: 21.6°C 
Se ee ee 
# (0 (mV) Data # (0 (mV) Data 
| 3 012) | 9S | 2 | (180) | 96 | 8 
[so | 467 | 009 | 16 G3 | 309 | 363 


Mode Shape for Length L/32, Area Ratio 0.7, and f=426.9 
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